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ABSTRACT 

We show that Supernova (SN) 2002kg in NGC 2403, initially classified as Type II- 
narrow (Iln), has photometric and spectroscopic properties unlike those of normal SNe. 
Its behavior, instead, is more typical of highly massive stars which experience the short- 
lived luminous blue variable (LBV) phase toward the end of their lives. The star, in 
fact, most resembles the LBV S Doradus in outburst. The precursor of SN 2002kg is 
the irregular, bright blue variable star 37 (V37), catalogued by Tammann & Sandage in 
1968. Using high-quality ground-based, multi-band images we can constrain the initial 
mass of V37 to be M ini ^ 40 M . We find that, although the spectra indicate a nitrogen 
enhancement, possibly revealing the products of CNO processing by V37 in the ejecta, 
the star lacks a substantial LBV nebula. The outburst from SN 2002kg/V37 is not 
nearly as energetic as the giant eruptions of the rj Carinae-like variables, such as SN 
1954J/V12, also in NGC 2403. SN 2002kg/V37, however, is among a growing number of 
"SN impostors" exhibiting a broad range of outburst energetics during a pre-SN phase 
of massive-star evolution. 

Subject headings: supernovae: general — supernovae: individual (SN 2002kg) — stars: 
early-type — stars: evolution — stars: variables: other — 
Hertzsprung-Russell (HR diagram) — galaxies: individual (NGC 2403) 



1. INTRODUCTION 

The evolution of the most massive stars is not 
well understood. After the main sequence, stars 
with masses M ^ 20-30 Mq should go through 
the red supergiant phase, or directly to a short- 
lived luminous blue variable (LBV) phase, to be- 
come Wolf-Rayet stars (e.g., Langer et al. 1994; 
Stothers & Chin 1996) before termination as su- 



pernovae (SNe; e.g., Woosley, Langer, & Weaver 
1993) or "collapsars" (e.g., MacFadyen & Woosley 
1999). LBVs which experience normal erup- 
tions or outbursts generally increase in bright- 
ness by one or two visual magnitudes, but re- 
main relatively constant in bolometric luminos- 
ity (Humphreys & Davidson 1994). However, the 
extremely rare case of rj Carinae (e.g., David- 



1 BASED IN PART ON OBSERVATIONS WITH THE NASA/ESA HUBBLE SPACE TELESCOPE, OBTAINED 
AT THE SPACE TELESCOPE SCIENCE INSTITUTE (STSCI), WHICH IS OPERATED BY AURA, INC., UNDER 
NASA CONTRACT NAS5-26555. 
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son & Humphreys 1997) shows that some very 
massive stars go through tremendous eruptive 
phases of mass ejection prior to the end of their 
lives. An estimate of the initial mass for rj Car 
is ^150 M Q (Hillier et al. 2001). The very mas- 
sive stars define the upper luminosity boundary of 
the Hertzsprung-Russell (HR) diagram. It would, 
of course, be very instructive to determine the 
prevalence of (and mass range for) the stars which 
undergo outbursts prior to becoming a SN. 

Several recent luminous events which have been 
identified as SNe are probably not genuine SNe 
at all. Instead, it is argued that these are more 
likely to be super-outbursts of very massive stars, 
analogous to r] Car (Goodrich et al. 1989; Filip- 
penko et al. 1995; Van Dyk et al. 2000). The 
energetics are comparable to those of SNe; hence, 
these objects act as SN "impostors." The "classi- 
cal" examples are SN 1954J/V12 in NGC 2403 
(Humphreys, Davidson, & Smith 1999; Smith, 
Humphreys, & Gehrz 2001; Van Dyk et al. 2005) 
and SN 1961V in NGC 1058 (e.g., Van Dyk, Fil- 
ippenko, & Li 2002, Van Dyk 2005; however, see 
Chu et al. 2004). Two more well-studied cases 
are SN 1997bs in M66 (Van Dyk et al. 2000) and 
SN 2000ch in NGC 3432 (Wagner et al. 2004). 
These objects were all underluminous and exhib- 
ited photometric behavior unlike that of any nor- 
mal SNe. For SN 1997bs, Van Dyk et al. (1999) 
identified a luminous precursor star, with My ~ 
—8.1 mag, in an archival Hubble Space Telescope 
(HST) WFPC2 F606W image. Unfortunately, 
no color information was available for the star 
to compare it to theoretical evolutionary models, 
in order to derive or strictly constrain the star's 
mass. 

In this paper we present another similar case, 
SN 2002kg, also in NGC 2403, which we simi- 
larly argue is not actually an explosion at the end 
of a massive star's life, but instead is the power- 
ful eruption of a LBV. SN 2002kg was discovered 
on 2003 October 22 (UT dates are used through- 
out this paper) by Schwartz & Li (2003) with 
the Tenegra II and Katzman Automatic Imaging 
Telescope (KAIT; Li et al. 2000; Filippenko et 
al. 2001; Filippenko 2005); it first became visible 
in images taken on 2002 October 26 (UT dates 
are used throughout this paper). 

In Figure 1 we show a KAIT image of the SN in 
the host galaxy. Filippenko & Chornock (2003) 
identified the SN spectroscopically as a Type II- 
"narrow" supernova (SN Iln; Schlegel 1990, Fil- 



ippenko 1991, 1997); it has Balmer emission lines 
exhibiting a narrow component atop a broader 
base, but with characteristics more similar to 
those of SN 1997bs (Van Dyk et al. 2000) than 
to the prototypical SNe Iln (e.g., SN 1988Z, Tu- 
ratto et al. 1993; SN 1995N, Fransson et al. 2002). 
In Van Dyk (2005) we first identified SN 2002kg 
with a previously known luminous irregular blue 
variable (V37) in the host galaxy and provided an 
estimate of the initial mass of this precursor star. 
Here we present a more detailed discussion. Weis 
& Bomans (2005) also have recently identified SN 
2002kg with V37 and discuss its historical light 
curve. 

2. OBSERVATIONS 
2.1. Photometry of the Outburst 
2.1.1. KAIT Monitoring 

Since the host galaxy has been monitored fre- 
quently by KAIT we are able to obtain unfiltered 
(~R) photometry for the SN 2002kg outburst over 
a period of more than one year. Li et al. (2003) 
discuss how the unfiltered magnitudes for an ob- 
ject can be transformed to standard R magni- 
tudes, with knowledge of its detailed color evo- 
lution. If the object's color remains unchanged, 
then the transformation is merely a constant off- 
set. We find that the SN has approximately the 
same color (within ~0.1 mag) at both an early 
and late time in the outburst (see §5); we could 
therefore follow the prescription in Li et al. (2003) 
to derive the offset. 

Equivalently, we can simply adjust the KAIT 
magnitudes, by about —0.7 mag, based on the 
offset between the unfiltered magnitude from the 
March 25.17 (JD 2452723.67) KAIT observation 
and the i?-band magnitude from the March 26.15 
(JD 2452724.65) Palomar observation (see §2.1.2; 
we consider the likely change in the R magni- 
tude over about one day to be insignificant for 
our purposes). We note that systematic errors in 
the magnitudes probably still exist for this KAIT 
light curve, due to the broad spectral range of the 
unfiltered images and the assumption of constant 
color for SN 2002kg (which may not apply for the 
entire outburst). We list the resulting R magni- 
tudes in Table 1 and show the light curve in Figure 
2, after further correcting for both extinction (see 
§3) and the true distance modulus to NGC 2403 
(fi = 27.48 mag; Freedman et al. 2001). 

The light curve shows highly erratic and un- 
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usual behavior for a SN. The maximum R lu- 
minosity is about —10.2 mag, more than 7 mag 
fainter than a typical SN Iln (see Van Dyk et 
al. 2000, their Figure 5). In fact, this light 
curve is more qualitatively similar to those of 
the possible n Car analog, SN 2000ch (Wagner 
et al. 2004, their Figure 3), or the LBV S Do- 
radus (Humphreys & Davidson 1994; their Figure 
2). Schwartz & Li (2003) note that SN 2002kg is 
not detectable, R ^ 19.5 mag (M R ^ —8.0 mag), 
in a KAIT image from 1998 November 13. Ap- 
parently we first caught the SN with KAIT while 
it was already undergoing the observed outburst 
(this is supported by the lack of detection at ear- 
lier epochs, based on data from other telescopes; 
see §2.1.2). The absolute magnitude of this out- 
burst was already ^2 mag brighter at that time 
than when the precursor star was in quiescence 
(see also §3). 

As Schwartz & Li (2003) point out, the light 
curve shows a brightening trend between 2002 Oc- 
tober 26, when the outburst was first noticed, and 
2003 early January. The light curve then shows a 
noticeable dip after 2003 January 1, before reach- 
ing maximum near 2003 mid-March. The bright- 
ness declines yet again by late March. No cover- 
age exists from KAIT between 2003 late March 
and mid-October, so we have no knowledge of the 
outburst behavior during that interval. However, 
by October 11 SN 2002kg has approximately the 
same brightness as it did in late March. During 
2003 October and November the object remains 
at a relatively constant brightness. We then see 
a pronounced dip in the light curve, after 2003 
late November through 2004 early March. The 
light curve rises again into 2004 early April, at 
which time our KAIT coverage terminates. This 
is similar to the behavior of SN 2000ch; Wagner 
et al. (2004) speculate that the dip seen for SN 
2000ch may be due to either a dust formation 
event, an eclipse, or changes within an optically 
thick ejected envelope. For SN 2002kg, this dip 
is much broader in time than was seen for SN 
2000ch. We suspect that for SN 2002kg the light 
curve dip is more likely due to opacity changes in 
the ejected envelope than to dust formation (see 
§4). 



2.1.2. Other Ground- Based Imaging 

As already mentioned, we obtained R, as well 
as BVI, images of SN 2002kg with the Palomar 
1.5-m telescope on 2003 March 26.15. In Figure 
3 we show the F-band image. We extracted in- 
strumental magnitudes for the relatively bright 
SN 2002kg in the Palomar images via fitting an 
empirically derived model point-spread function 
(PSF) for each band, using DAOPHOT (Stetson 
1987) within IRAF 2 . Unfortunately, no photomet- 
ric calibration was obtained during these Palo- 
mar observations. However, Larsen & Richtler 
(1999) also obtained UBVRI (plus Ha) images 
of NGC 2403 with the Nordic Optical Telescope 
(NOT) 3 on 1997 October 13, prior to SN 2002kg, 
under very good seeing conditions of 0'.'8. (We will 
discuss these NOT images in more detail in §3.) 
We extracted instrumental magnitudes from the 
NOT broad-band images using IRAF, first with a 
4" aperture and then with a PSF for each band. 
Calibration was established using photometry of 
a number of isolated stars through a 20" aperture, 
which matches the aperture used to establish the 
calibration through standard-star observations by 
Larsen (1999; we also applied the appropriate ex- 
tinction corrections, supplied via a private com- 
munication from S. S. Larsen). We have com- 
pared our photometry with the BV bright-star 
photometry of Sandage (1984) and UBVR pho- 
tometry of Zickgraf & Humphreys (1991) and find 
good agreement across all bands. 

We tie the photometry from the Palomar im- 
ages to that from the calibrated NOT images and 
list the results in Table 2. As mentioned above in 
§2.1.1, we used the R magnitude for SN 2002kg 
from these Palomar data to adjust the unfiltered 
KAIT light curve. We also include this R mag- 
nitude from Palomar in the light curve shown in 
Figure 2. 

A number of images of NGC 2403 obtained 
in several bands prior to 2002 October are also 
available in the Isaac Newton Group archive. 
These images have a range of sensitivities, and 
SN 2002kg is not detected in any of them. We 
list the resulting upper limits to detection of the 
outburst in Table 2. 



2 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, which 
are operated by the Association of Universities for Research in Astronomy Inc., under cooperative agreement with the 
National Science Foundation. 

3 The NOT images were contributed to the NASA/IPAC Extragalactic Database (NED) and posted for public distri- 
bution on their website, http://nedwww.ipac.caltech.edu. 
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2.1.3. HST Imaging 

Additionally, SN 2002kg was observed on 2004 
August 17 with the Wide Field Channel (WFC) 
and on 2004 September 21 with the High Res- 
olution Channel (HRC) of the Advanced Cam- 
era for Surveys (ACS), on-board HST. The Au- 
gust observations were obtained during a cam- 
paign by program GO-10182 (PI: Filippenko) to 
observe SN 2004dj, the third historical SN in NGC 
2403. The imaging strategy for SN 2004dj was 
purposely designed so that the site of SN 2002kg 
would also be contained in the WFC images. The 
September observations were obtained as part of 
a larger Snapshot program observing SNe (GO- 
10272; PI: Filippenko). The WFC exposures were 
600 s in F475W, 650 s in F658N, and 350 s each 
in F606W and F814W. The HRC exposures were 
840 s in F435W and 360 s in F625W. The im- 
ages were all extracted from the HST archive. 
For the WFC data we recombined the individual 
"fit" exposures in all four bands at the Space Tele- 
scope Science Institute (STScI), in order to more 
reliably remove cosmic-ray hit features from the 
images, and, in effect, produce our own drizzled 
(Fruchter & Hook 2002) "drz" images indepen- 
dent of the HST data pipeline. 

The WFC F606W image is shown in Figure 
4. We performed PSF-fitting photometry within 
IRAF/DAOPHOT on the broad-band WFC im- 
ages. The empirical PSF for each band was de- 
rived from the ~ 10 brightest, uncrowded, and 
unsaturated stars in these images, based on a 
0'.'5 aperture. For the WFC narrow-band image, 
we performed aperture photometry only, using a 
0'.'5 aperture. We determined that any correction 
to the magnitudes for the degradation of charge- 
transfer efficiency (CTE) is <0.2% in all bands 
(see Riess & Mack 2004). Next, for the F475W 
and F606W filters we applied the corrections from 
a (X'5 aperture to infinite aperture from Table 5 
in Sirianni et al. (2005); for the F814W filter, we 
applied the correction from their Table 6, assum- 
ing the effective wavelength from their Table 8 for 
approximate spectral type F2. We then applied 
the zero points (VEGAMAG for the broad bands, 
and STMAG for the narrow band) in Table 10 of 
Sirianni et al. to the corrected instrumental mag- 
nitudes. The resulting flight system magnitudes 
are m F475W = 19.43±0.02, m F606W = 19.26±0.02, 
"T-F814W — 19.13±0.03, and mp658N = 18.67±0.14. 

4 IDL is the Interactive Data Language, a product of R 



In the HRC images SN 2002kg is relatively 
isolated, so we simply employed aperture pho- 
tometry (assuming a / .'5 radius aperture) within 
IRAF. We then applied the corrections to infi- 
nite aperture from Table 5 and the zero points 
from Table 11 in Sirianni et al. (2005). The re- 
sulting flight system magnitudes are mp435w = 
19.41 ± 0.05 and m F625 w = 19.16 ± 0.04. 

Finally, we transformed the broad-band magni- 
tudes into the Johnson-Cousins system, also fol- 
lowing Sirianni et al., keeping in mind the fact 
that SN 2002kg shows strong emission lines (Fil- 
ippenko & Chornock 2003, and also §2.2), whereas 
the stars used to establish the transformations all 
show normal absorption-line spectra. We list the 
resulting magnitudes in the standard photomet- 
ric system in Table 2. We also include the HRC 
i?-band data in the light curve shown in Figure 2. 

2.2. Spectroscopy of the Outburst 

Optical spectra of SN 2002kg were obtained at 
the Keck 10-m telescope with the Low Resolution 
Imaging Spectrometer (LRIS; Oke et al. 1995) 
on 2003 January 6 (the confirmation spectrum; 
JD 2452646.0) during the early brightening trend, 
on 2003 November 29 (JD 2452973.0), prior to 
the pronounced dip in brightness, and on 2004 
December 12 (JD 2453352.03) at relatively later 
times. These observations had spectral resolution 
of ~6 A across most of the spectral range, 3000- 
9430 A (the resolution is somewhat worse in the 
near- infrared) . The 2004 December spectrum was 
obtained with a 600-line grism on the blue side, 
whereas the other two LRIS spectra were obtained 
with a 400-line grism in the blue. At Ha the reso- 
lution (FWHM) of these spectra is ~320 km s -1 . 

In addition, a spectrum was obtained on 2003 
November 20 (2452964.0), while the object was 
still bright, with the blue channel of the spec- 
trograph on the 6.5-m Multiple-Mirror Telescope 
(MMT). A 300-line grating was used, with spec- 
tral range 3700-8300 A. For the Keck spectra con- 
ventional data reductions were performed within 
IRAF, including bias subtraction, flat-fielding, 
and wavelength calibration, while our own IDL 4 
routines were used for flux calibration and re- 
moval of telluric lines (e.g., Matheson et al. 2000). 
Reduction of the MMT spectrum was performed 
in a similar fashion. 

The spectra are shown in Figure 5. The over- 
all shape of the continuum has remained rela- 

ch Systems, Inc. 
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tively blue throughout 2003 and 2004, peaking 
near 3700 A. The spectra all show several strong 
emission lines, primarily of H, but also of possible 
weak Fe II at 4500-4600 A and at 5150-5350 A, 
as well as possible blends of H with He I. Figure 
6 shows just the Ha plus [N II] AA6548, 6584 A 
lines. The narrow Balmer profiles are unresolved, 
and the broad profiles at Ha remain relatively 
constant throughout 2003 and 2004, with Gaus- 
sian velocity widths a m 370 km s _1 . Such a 
velocity width for the line is uncharacteristically 
low for a young SN, but consistent with the ter- 
minal wind velocities for LBVs. The presence of 
the prominent, unresolved [N II] emission lines, 
while other forbidden lines (e.g., due to O and 
S) appear weak or absent, suggests a possible N 
enhancement in the circumstellar gas. 

From the line profiles in Figure 6 it is clear that 
the power in the Ha+[N II] emission has steadily 
declined. This trend is consistent with that of the 
-R-band light curve for the observation dates of the 
spectra. In fact, integrating over the continuum- 
subtracted Ha+[N II] line profiles shown in Figure 
6, we find an observed flux Fh q+ [nii] ~ 1.3 x 10 _ 14 
erg cm -2 s _1 (uncorrected for extinction) in 2003 
January, ~ 1.1 x 10~ 14 erg cm' 2 s -1 in 2003 
November, and ~ 4.5 x 10 -15 erg cm -2 s -1 in 
2004 December. (The absolute flux scale for the 
Keck spectra is accurate to ~ 20% or better.) 

In the blue portion of the spectra, shown in 
Figure 7, we see most notably the presence of 
various absorption lines, mostly due to higher- 
order Balmer lines, and possible blends with He I. 
Strong Ca II H and K lines (the H component is 
also probably blended with He) are also seen, par- 
ticularly in 2003, but we suspect these lines are 
dominated by interstellar absorption local to the 
SN (the extinction toward the SN appears to be 
significant; see §3). 

Such absorption features are not typically seen 
in SN spectra. (We have been very careful in the 
extraction of the spectra from the 2-dimensional 
images and can safely rule out contamination 
by a neighboring star.) The equivalent widths 
of these lines indicate an underlying star with 
spectral type and luminosity class of a late O 
or early-to-mid B supergiant. We compare the 
SN 2002kg spectra to those of the LBV S Dor 
(from 1996, during outburst; Massey 2000) and 
the well-studied S Dor-type LBV AG Carinae (in 
a deep minimum state; Walborn & Fitzpatrick 
2000). Although the Balmer emission lines for 



SN 2002kg lack the P-Cygni profiles seen in the 
AG Car spectrum, the SN spectra, overall, show a 
remarkable resemblance to that for S Dor in out- 
burst. We therefore consider the spectroscopic 
evidence quite compelling that SN 2002kg is not 
a true SN, but instead a S Dor-like LBV in out- 
burst. A similar conclusion was also reached by 
Weis & Bomans (2005). 

3. VARIABLE 37: THE VERY MASSIVE 
PRECURSOR OF SN 2002KG 

From the Palomar images (§2.1.2) we mea- 
sure an accurate absolute position for the SN of 
a(J2000) = 7 h 37 m 01166, 5(J2000) = 65° 34' 28"0, 
with a total uncertainty of ±0 / .'2, using the Two 
Micron All Sky Survey (2MASS) as the astro- 
metric reference (see Van Dyk, Li, & Filippenko 
2003). This position differs by l'!7 from that re- 
ported by Schwartz & Li (2003), and results in an 
offset from the host galaxy's nucleus of 63'.'6 east 
and 10172 south. We then applied the same astro- 
metric grid to the NOT images (with a ±0?2 total 
uncertainty in the grid). In Figure 8 we show the 
NOT K-band image from 1997. When the abso- 
lute SN position is superimposed on this image, 
it coincides exactly with the star indicated by the 
arrow in the figure. This star is almost certainly 
the precursor of SN 2002kg. 

We find that this star had V = 20.63 ± 0.04, 
U - B = -0.96 ± 0.05, B -V = -0.08 ± 0.05, 
V - R = 0.09 ± 0.04, and V - I = -0.06 ± 0.05 
mag (see Table 2). (Uncertainties include the 
uncertainties in the standard calibration in each 
band from Larsen 1999, his Table 1, added in 
quadrature with the measurement uncertainties 
from IRAF.) The fainter (V = 21.01 ±0.04 mag), 
redder (V - I = 0.50 ± 0.05 mag) star ~ l'.'O 
to the northeast is well outside the total uncer- 
tainty in the SN position in the NOT image, 073 
(summing in quadrature the total uncertainty in 
the SN position measured from the Lick image 
and the uncertainty in the astrometric grid on 
the NOT image). Also, this fainter star appears 
at approximately this same brightness (to within 
one magnitude) and color in the HST images as in 
the NOT images, whereas SN 2002kg (and there- 
fore its precursor) have changed more significantly 
both in brightness and in color (compare Figure 8 
with Figure 3 and see Table 2). We note that Weis 
& Bomans (2004) also attempt to identify the 
precursor, but, using pre-outburst ground-based 
images with inferior spatial resolution, they are 
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unable to distinguish the star we have identified 
as the precursor from its fainter neighbor without 
consulting the post-outburst HST images. 

This precursor star had actually been identified 
first by Tammann & Sandage (1968) in their pho- 
tographic study of supergiants in NGC 2403. We 
have determined this from a careful astrometric 
comparison of their Figure 5, an enlargement of a 
103aO photographic plate of the galaxy, with the 
NOT V-band image. The SN precursor star we 
have identified lines up very well with the position 
of the catalogued bright blue irregular variable 
star, no. 37 (hereafter, V37) from Tammann & 
Sandage (1968). We consider the positional coin- 
cidence a compelling indication that the two stars 
are one and the same (see also Weis & Bomans 
2005). 

Tammann & Sandage (1968) describe the B 
light curve for V37 as exhibiting relatively small 
amplitudes of ~ 2 mag, between about 19 and 21 
mag, and characteristic fluctuation timescales of 
several years (see their Figure 10). Where data 
were available between about 1948 and 1963, the 
observed B — V color of V37 apparently fluctu- 
ated between +0.20 and —0.36 mag. The maxi- 
mum absolute B magnitude for V37 during their 
monitoring was ~ —8.8, although generally it was 
~ —6.5 mag. 

We can estimate the extinction and redden- 
ing to this star by placing it in color-color dia- 
grams (Figure 9), based on the NOT and HST 
ACS/WFC photometry. Also shown in both di- 
agrams are the loci of the main sequence and 
supergiants (Drilling & Landolt 2000), and the 
reddening vector for OB supergiants, assuming 
a Cardelli, Clayton, & Mathis (1989) reddening 
law. Figure 9a, in particular, implies that the star 
is a blue supergiant reddened along this vector, at 
A v = 0.6±0.2 mag (E[B-V] = 0.19±0.06 mag). 
Most of this extinction must be in the star's local 
environment; the component of the Galactic fore- 
ground extinction alone is only Ay = 0.13 mag 
(Schlegel, Finkbeiner, & Davis 1998). In Figure 9 
we also show the loci of stars in the precursor's im- 
mediate (~ 200 pc) environment; the reddening 
in the environment appears to be quite variable, 
with many stars more reddened than the precur- 
sor. (The fainter star ~ 1" to the northeast is 
quite reddened, with Ay ^ 2 mag.) The extinc- 
tion toward the precursor is generally consistent 
with the interstellar extinction for stars in its en- 
vironment. 



Correcting the star's observed photometry for 
the implied extinction and reddening, as well 
as for the true distance modulus of NGC 2403 
(/i = 27.48 ±0.1 mag; Freedman et al. 2001), we 
find My = -7.45±0.23, (U-B) = -1.02±0.07, 
(B-V) = -0.13±0.08, (V-R) = -0.01±0.05, 
and (V - I) = -0.30 ± 0.09 mag. (The total un- 
certainty in the absolute magnitude is the mea- 
surement uncertainty and the uncertainty in the 
distance modulus, added in quadrature; the to- 
tal uncertainties in the colors include the mea- 
surement uncertainties and the uncertainty in the 
reddening, assuming the Cardelli et al. 1989 red- 
dening law, again added in quadrature.) The pre- 
cursor indeed appears to be a very luminous, blue 
supergiant of type 09 to B2. This identifica- 
tion of the star's type, based on color alone, is 
consistent with that inferred from the underlying 
absorption-line features in the SN spectra. 

During outburst, as seen in Figure 9b from the 
2003 March Palomar data, the star has developed 
the observed colors of a late A-type or early F- 
type supergiant. V37, therefore, is following the 
typical behavior of a normal LBV during out- 
burst: a "pseudo-photosphere" expands in an op- 
tically thick stellar wind, increasing the star's V 
brightness by sim2-3 mag (from V = 20.63 to 
V ~ 18), while decreasing the apparent stellar 
temperature to T eS « 8000 K (although the con- 
cept of T c fj, as applied to LBVs in eruption, loses 
its true meaning, and really something more like 
an "apparent temperature" T app is appropriate; 
see Humphreys & Davidson 1994), without a sub- 
stantial increase in reddening due to dust forma- 
tion or change in the bolometric magnitude, M^ Q \ 
(Humphreys & Davidson 1994). The formation 
of extended, dense, dusty ejecta, such as in the 
super-outbursts of 77 Car (Humphreys & David- 
son 1994; Davidson & Humphreys 1997) and SN 
1954 J/ V12 (Van Dyk et al. 2005), appears to be 
relatively exceptional behavior. 

At T eff = 7000-9000 K, the bolometric correc- 
tion is zero or nearly zero (at most, about —0.2 
mag). Therefore, for V37 at maximum light dur- 
ing the outburst, Mboi ~ My —9.8 mag (as- 
suming Ay = 0.6 mag and that the 2003 March 
Palomar V measurement is near maximum, to 
within ~ 0.2 mag). If we assume that Mboi has re- 
mained constant, pre-outburst and at maximum, 
then, with the intrinsic pre-outburst color of the 
star, (B — V) = —0.13 mag, and the correspond- 
ing T eS S3 15500 K, we can place the star in a 
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HR diagram and estimate the initial mass of the 
precursor. 

We show the star's location in such a HR dia- 
gram in Figure 10. For comparison we show the 
theoretical evolutionary tracks from Lejeune & 
Schaerer (2001), with enhanced mass loss for the 
most massive stars, for a 120, 85, 60, and 40 Mq 
star. In Figure 10a we show the tracks assuming 
solar metallicity (Z = 0.02). It is also possible, 
however, that the star may have had abundances 
different from solar. SN 2002kg is 11$! 5 from the 
galaxy nucleus, which translates to a galactocen- 
tric distance R xs 1.9 kpc. From a study of the 
metallicity gradient in NGC 2403 by Garnett et 
al. (1997), the value of log(0/H) + 12 at this dis- 
tance is 8.6 dex. With the solar log(0/H)+12 
abundance at 8.8 dex (Grevesse & Sauval 1998), 
the metallicity in this region of the galaxy may ac- 
tually be somewhat subsolar; thus, in Figure 10b 
we show the tracks calculated at the next lowest 
metallicity, Z = 0.008, as well. 

What can immediately be seen from this figure, 
regardless of metallicity, is that the SN 2002kg 
precursor, V37, is well off the zero-age main se- 
quence and is in a more evolved evolutionary 
phase. We conclude that the initial mass for V37 
is Mzams ^ 40 Mq. Unfortunately, we cannot 
be more restrictive in our initial mass estimate, 
since, due to the uncertainties in the photometry, 
distance, and reddening of the star, the termini 
of the 120, 85, and 60 Mq tracks are also consis- 
tent with this photometry. Nonetheless, the HR 
diagrams indicate that V37 is indeed a very mas- 
sive star in the upper regions of the HR diagram, 
consistent with the theoretical mass estimates of 
stars which go through the LBV phase. 

4. A LBV NEBULA AROUND SN 2002KG/V37? 

From the NOT and ACS/WFC Ha images, we 
can investigate the gaseous environment of V37, 
both before and during outburst. Weis & Bomans 
(2005) have pointed out that in both the HST 
images obtained in the F658N band and ground- 
based images from 2001, net Ha emission is seen 
at the position of the star. However, Humphreys 
& Aaronson (1987) describe a spectrum of the 
star obtained in 1985 as consisting of absorption 
lines and a blue continuum (consistent with fea- 
tures seen in our own spectra of SN 2002kg/V37), 
but with no emission lines visible. Weis & Bo- 
mans speculate that the observed Ha emission is 
due to the formation of a LBV nebula, associated 



with the current outburst or possibly ejected dur- 
ing an earlier evolutionary phase (e.g., Lamers et 
al. 2001). 

We confirm that net Ha (+ [N II]) emission is 
associated with the star at least as early as 1997, 
by subtracting the NOT R continuum image from 
the Ha narrow-band image, after first intensity- 
scaling the two images (this technique can do a 
reasonable job of removing the continuum, when 
a specific narrow-line bandpass filter contiguous 
to Ha in wavelength is not available; see Waller 
1990). We show the result of this subtraction in 
Figure 11; there is a very compact (i.e., unresolved 
and not extended) Ha (+ [N II]) source at the ex- 
act position of the star. 

We would expect a star with spectral type 09- 
B2 to have an associated photoionization (H II) 
region. If we scale the NOT net Ha (+ [N II]) im- 
age by fluxes for the giant H II regions S158, S256, 
and S298 in NGC 2403 (Sivan et al. 1990; Drissen 
et al. 1999), we find a flux F Hq+ [nii] ~ 1.1 x 10 -15 
erg cm~ 2 s _1 for the nebula (uncorrected for ex- 
tinction). We can correct for extinction using 
the relation A~n a = 0.81^4y from Viallefond & 
Goss (1986), and with Ay = 0.6 mag (assuming 
the line-of-sight extinction to the H II region is 
the same as that toward the precursor), we find 
1.7 x 10~ 15 erg cm" 2 s" 1 . At the distance of NGC 
2403 (Freedman et al. 2001) this results in a lumi- 
nosity Lhch-[nii] = 2.1 x 10 36 erg s -1 . Using the re- 
lation A^L yc = 7.25 x 10 n LH a to derive the number 
of ionizing Lyman continuum photons (1.5 x 10 48 ), 
and using Table II of Panagia (1973), we find that 
this is the approximate equivalent of the ionizing 
flux from a single B0 I star. That is, the Ha emis- 
sion seen in 1997 is at least consistent with a H II 
region produced by the ionizing flux of a single 
massive blue supergiant — the outburst precur- 
sor, V37. 

Subtracting the continuum in the F606W band 
from the 2004 August ACS/WFC F658N image, 
and scaling the net emission image to the fluxes 
of the giant H II regions S256 and S298, results in 
a flux F Hq+ [nii] ~ 5.2 x 10~ 15 erg cm -2 s -1 (un- 
corrected for extinction) for SN 2002kg. This is 
approximately a five-fold increase in the Ha flux 
during outburst in late 2004, relative to the pre- 
outburst state. We have already seen from §2.2 
that in 2003 January, prior to maximum light, 
the star was already in outburst and was emit- 
ting ~ 12 times more Ha flux than in its qui- 
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escent state. This 2004 August flux from the 
HST image is consistent with the declining fluxes 
derived from the integrated spectral line profiles 
from 2003 through 2004 (see §2.2). 

The outburst Ha flux is likely originating in the 
dense, optically thick wind from the star. But, 
can we really consider this a LBV nebula? The 
relatively strong [N II] lines do imply the pres- 
ence of circumstellar g does the indication 
of CNO-processed material in the wind. However, 
a normal LBV eruption does not generally result 
in the ejection of a massive shell, or nebula, in 
the first place. If SN 2002kg is similar to S Dor, 
then we might expect any nebula to be quite com- 
pact. Based on HST imaging, Weis (2003) places 
an upper limit on the S Dor nebula diameter of 
0.25 pc. Similarly, we find that the source is un- 
resolved (£1.8 WFC1 pixels FWHM) in the 2004 
August HST F658N image which, at the distance 
to NGC 2403 (3.13 Mpc; Freedman et al. 2001), 
implies a diameter £1.4 pc for the nebula. If the 
wind has been expanding at a constant 370 km 
s- 1 between 2002 Oct 26 and 2004 Aug 17, this 
implies that the wind-blown shell, or nebula, is 
only ~ 0.001 pc in diameter. We infer that the 
detected Ha emission is from the ejecta of the 
recent eruption. 

5. DISCUSSION AND CONCLUSIONS 

We have shown that SN 2002kg in NGC 2403 is 
unlike a typical Type II, core-collapse SN, given 
both its spectroscopic and photometric proper- 
ties. Although spectroscopically SN 2002kg re- 
sembles the SNe Iln, the ejecta expansion velocity, 
~370 km s" 1 , is much lower than what is observed 
for the prototypical SNe Iln (e.g., SN 1988Z; Tu- 
ratto et al. 1993). Additionally, the star increased 
only ~ 2 mag in V. We deduce, therefore, that 
SN 2002kg is, in fact, not a genuine SN (strictly 
defined as the explosion of a star at the end of its 
life), but instead represents an eruption of a nor- 
mal, or "classical," supergiant LBV star, possibly 
analogous to S Dor. With Mboi ~ — 9.8 mag and 
logT eff ps 4.19 for the precursor in its quiescent 
state, the star clearly lies on the "S Dor insta- 
bility strip" (Wolf 1989; Humphreys & Davidson 
1994). 

We further support this conclusion by show- 
ing in Figure 12 the absolute light curve for SN 
2002kg, relative to the light curves for other SN 
impostors (SNe 1961V, 1954J/V12 1997bs and 
2000ch) in a similar vein as Van Dyk et al. (2000; 



their Figure 6). Although the bandpasses through 
which these light curves were generated differ 
from one object to the other, the relative similar- 
ities in luminosity, energy, and variability argue 
strongly that SN 2002kg is also a SN impostor. 

SN 2002kg appears from its light curve to be 
the least energetic example of a SN impostor 
known so far. SNe 1954J/V12, 1961V, and rj Car 
may represent the high-energy end of a pre-SN 
super-outburst phase, with SN 1997bs and SN 
2000ch somewhere in between. The lower ener- 
getics can also be seen in a comparison (Fig. 13) 
of the early-time spectra of SNe 2002kg, 1997bs, 
and 2000ch. The overall strength of the Balmer- 
line profiles in SN 2002kg is lower than in the 
other two SNe. Moreover, the bases of the pro- 
files in SN 2002kg are not nearly as broad as 
in SNe 1997bs and 2000ch, consistent with the 
lower ejecta velocities in SN 2002kg. (SN 2000ch 
also shows more of a true P-Cygni profile, with 
shallow blue absorption troughs in the Balmer 
lines.) SN 2002kg is clearly less energetic than SN 
1954 J/ V12 in the same host galaxy, NGC 2403; 
this is reflected both in the lower expansion ve- 
locity (~ 370 km s^ 1 , versus ~ 700 km s _1 at age 
~50 yr; Van Dyk et al. 2005), and fainter appar- 
ent brightness (my [max] ps 18.3 mag versus 16.1 
mag; Humphreys & Davidson 1994). 

The precursor of SN 2002kg is the irregular 
blue variable star V37 in NGC 2403, identified by 
Humphreys & Aaronson (1987) and Humphreys 
& Davidson (1994) as a LBV. We have also di- 
rectly placed constraints on the initial mass of 
this precursor to be ^ 40 Mq. As a result of the 
outburst, the star became redder and attained a 
nearly constant color (the B — V color remained 
relatively unchanged for nearly 1.5 yr; the over- 
all constant shape of the spectral continua fur- 
ther supports this), as the pseudo-photosphere of 
the optically thick envelope expands. The abso- 
lute flux in the Ha line decreased between early 
2003 and late 2004; this is also reflected in the 
B — R color of the star (0.41 mag in 2003 March 
and 0.26 mag in 2004 August /September). The 
constant color argues against any dust formation 
event for SN 2002kg/V37, such as is evident in the 
Homunculus around r\ Car and in SN 1954 J/ VI 2 
(Van Dyk et al. 2005), although dust might form 
later for the SN 2002kg/V37 ejecta. 

SN 2002kg/V37 is one of a growing number of 
objects which "impersonate" SNe Iln. These im- 
postors exhibit a broad range of properties and 
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actually represent a pre-SN evolutionary phase of 
very massive stars. It is imperative to continue 
monitoring the star to detail its further evolu- 
tion. Other recent, less well-studied cases, such 
as SN 1999bw in NGC 3198, SN 2001ac in NGC 
3504, and SN 2003gm in NGC 5334, also require 
detailed study, to provide valuable additional ex- 
amples of this rare phenomenon and improve our 
understanding of massive-star evolution. 
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Table 1 

KAIT Photometry of SN 2002kg 



u i elate 


J-D — Z40UUUU 


K (mag) 


1998 Nov 13 


1130 


>19.5 


2002 Oct 26.56 


2574.06 


18.4(0.1) 


2002 Nov 01.53 


2580.03 


18.3(0.1) 


2002 Nov 12.54 


2591.04 


18.1(0.2) 


2002 Nov 16.44 


2594.94 


18.4(0.1) 


2002 Nov 20.40 


2598.90 


18.4(0.2) 


2002 Dec 02.49 


2610.99 


18.0(0.1) 


2003 Jan 08.40 


2647.90 


17.9(0.1) 


2003 Jan 16.29 


2655.79 


18.0(0.1) 


2003 Jan 26.34 


2665.84 


18.2(0.1) 


2003 Jan 31.27 


2670.77 


18.2(0.1) 


2003 Feb 08.26 


2678.76 


18.1(0.1) 


2003 Feb 23.23 


2693.73 


18.1(0.1) 


2003 Mar 03.20 


2701.70 


18.0(0.1) 


2003 Mar 08.17 


2706.67 


18.1(0.1) 


2003 Mar 13.21 


2711.71 


17.8(0.1) 


2003 Mar 25.17 


2723.67 


18.0(0.1) 


2003 Oct 11.49 


2923.99 


18.1(0.1) 


2003 Oct 18.54 


2931.04 


18.3(0.1) 


2003 Oct 26.54 


2939.04 


18.0(0.1) 


2003 Nov 10.53 


2954.03 


18.1(0.1) 


2003 Nov 23.48 


2966.98 


18.1(0.1) 


2004 Feb 01.32 


3036.82 


18.5(0.1) 


2004 Feb 11.26 


3046.76 


18.8(0.1) 


2004 Mar 03.19 


3068.69 


19.1(0.3) 


2004 Mar 08.13 


3073.63 


19.4(0.2) 


2004 Mar 17.16 


3082.66 


18.7(0.1) 


2004 Apr 04.20 


3099.70 


18.4(0.1) 



Uncertainties in the magnitudes are given in parentheses. 
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Table 2 

SN 2002kg Photometry from Other Sources 



UT date 


U 


B 


V 


R 


I 


Source a 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




1987 Feb 04 




£20.3 


£20.4 


£20.1 




JKT 


1992 Mar 22 




£16.0 


£15.8 


£15.7 




INT 


1997 Oct 13 


19.59(04) 


20.55(03) 


20.63(04) 


20.54(02) 


20.69(03) 


NOT 


1998 Oct 05 




£20.5 




£20.4 




JKT 


1999 Apr 25 




£20.1 




£20.1 




JKT 


2003 Mar 26 




18.45(03) 


18.32(04) 


18.04(02) 


17.97(02) 


P60 


2004 Aug 17 b 




19.58(05) 


19.39(06) 




19.21(06) 


HST 


2004 Sep 21 b 




19.41(07) 




19.13(04) 




HST 



Note: uncertainties of hundredths of a magnitude are indicated in parentheses. 

a JKT=Johannes Kepler Telescope; INT=Isaac Newton Telescope; NOT=Nordic Optical 
Telescope; P60=Palomar 1.5-m telescope; HST=Hubble Space Telescope. 

b Transformed from the flight system magnitudes via synthetic photometry of normal stars 
with SYNPHOT; see text. 



Fig. 1. — SN 2002kg (indicated by tick marks) in NGC 2403 in an unfiltered (~ R) image obtained with 
the Katzman Automatic Imaging Telescope (KAIT) 0.76-m telescope on 2002 October 26. The image 
is Q'.Q x Q'.Q in size. North is up, and east is to the left. 
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Fig. 2. — The unfiltered (~ i?-band) light curve {open squares) of SN 2002kg obtained using KAIT. 
(See Table 1; see text regarding the photometric zero point and systematic errors in the light curve.) 
Also shown are the i?-band data for the SN obtained with the Palomar 1.5-m telescope on 2003 March 
26.15 {filled square) and with the HST ACS on 2003 September 21.07 {filled triangle). The photometry 
has been corrected for extinction (see text) and for the true distance modulus of the host galaxy, NGC 
2403 (Freedman et al. 2001). 
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Fig. 3. — SN 2002kg (indicated by the arrow) in a ^/-band image obtained at the Palomar 1.5-m 
telescope on 2003 March 26.15. 



Fig. 4. — SN 2002kg (indicated by the arrow) in a F606W image obtained with the ACS/WFC on-board 
HST on 2004 August 17.11. This figure is to approximately the same scale and orientation as Figure 3. 
Note that the SN field is near the edge of the ACS image. 
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Fig. 5.— Spectra of SN 2002kg obtained at the 10-m Keck telescope (2003 Jan 6, 2003 Nov 29, and 
2004 Dec 12) and the 6.5-m MMT (2003 Nov 20). 
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Fig. 6. — A comparison of the spectra shown in Figure 5 in the wavelength region containing Ha and 
the [N II] lines. The narrow emission component of Ha is unresolved, and the broad component remains 
relatively constant in width throughout 2003 and 2004. The prominent, unresolved [N II] AA6548, 6584 
A emission lines suggest the presence of N-enhanced circumstellar material. 
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3800 4000 4200 4400 4600 4800 

Rest Wavelength (A) 

Fig. 7. — Expanded views of the spectra of SN 2002kg shown in Figure 5 (redshift-corrected assuming 
z = 0.000437), in the blue wavelength region. Also shown for comparison are spectra of the LBVs AG 
Car (Walborn & Fitzpatrick 2000) and S Dor (Massey 2000; redshift-corrected assuming z = 0.000927). 
Various lines are labelled. 



Fig. 8. — The precursor of SN 2002kg (indicated by the arrow), as seen in a V'-band image of the host 
galaxy, NGC 2403, obtained with the 2.6-m Nordic Optical Telescope (NOT) on 1997 Oct 13 (Larsen 
& Richtler 1999). This figure is to nearly the same scale and orientation as in Figure 3. The location 
of the precursor is coincident with the irregular blue variable V37, reported by Tammann & Sandage 
(1968). 
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Fig. 9.— The (a) (U — B, B — V) and (b) (B — V, V — I) color-color diagrams of SN 2002kg and its 
stellar environment (within ~ 500 pc). The colors of the precursor star (filled circle) measured from 
the NOT images are shown in both diagrams. The open squares represent the colors of the stars in the 
environment measured from the NOT images for panel (a) and the HST ACS/WFC images for panel (b) 
(after transformation from flight system to standard Johnson-Cousins colors). Also shown are the loci 
for supergiant (long- dashed line) and main sequence (short-dashed line) stars, as well as the reddening 
vector (solid line) from Ay = to Ay = 3 mag, assuming the Cardelli et al. (1989) reddening law. From 
these diagrams, we infer Ay = 0.6 ± 0.2 mag for the precursor star. In (b) we show the colors of the SN 
during outburst from the 2003 Palomar data (filled square) and from the 2004 HST/ ACS data (filled 
triangle; the open triangle represents the SN colors after an assumed correction for contamination by 
the strong Ha emission in the F606W band; a dotted line connects the two symbols). 
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Fig. 9. — (Continued.) 
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Fig. 10. — HR diagrams, showing the locus for the precursor (five-pointed star), after correction for 
reddening and the true distance modulus of the host galaxy, NGC 2403, and assuming that Mt, i has 
remained constant before outburst and at maximum. Also shown for comparison are the model stellar 
evolutionary tracks from Lejeune & Schaerer (2001), with enhanced mass loss, for 120 Mq (dotted 
line), 85 Mq (short-dashed line), 60 Mq (solid line), and 40 Mq (long-dashed line). Panel (a) shows 
the tracks for solar metallicity (Z = 0.02), and panel (b) shows the tracks for sub-solar metallicity 
(Z = 0.008). 
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Fig. 10. — (Continued.) 
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Fig. 11. — The continuum-subtracted Ha image of the host galaxy, NGC 2403, obtained with the 2.6- 
m Nordic Optical Telescope (NOT) on 1997 Oct 13 (Larsen & Richtler 1999), showing net emission 
surrounding the precursor of SN 2002kg (indicated by the arrow). This figure is shown to the same 
scale and orientation as Figure 8. The nebula in this image is spatially unresolved, with diameter ^ 14 
pc (FWHM). The Ha luminosity is consistent with photoionization by a single blue supergiant star. 
Artifacts, due to incomplete continuum subtraction, are also seen in the figure. 
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Fig. 12. — A comparison of the absolute light curve for SN 2002kg in NGC 2403 (filled triangles and 
solid line; also see Figure 2) with that of rj Car (dotted line), SN 1954J in NGC 2403 (long-dashed, 
dotted line), SN 1961V in NGC 1058 (short- dashed, dotted line; adapted from Humphreys, Davidson, & 
Smith 1999), SN 1997bs in M66 (short-dashed line; Van Dyk et al. 2000), and SN 2000ch in NGC 3432 
(long-dashed line; Wagner et al. 2004). We note that the bandpasses through which each of these light 
curves were generated differed from one object to the other. 
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Fig. 13. — A comparison of the early-time optical spectrum of SN 2002kg (see Figure 5; redshift- 
corrected assuming z = 0.000437) with those of SN 1997bs in M66 (Van Dyk et al. 2000; redshift- 
corrected assuming z = 0.002425) and SN 2000ch in NGC 3432 (Wagner et al. 2004; redshift-corrected 
assuming z = 0.002055). 
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